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Abstract 

We develop an improved Statistical Multifragmentation Model that provides the capability to 
calculate calorimetric and isotopic observables with precision. With this new model we examine 
the influence of nuclear isospin on the fragment elemental and isotopic distributions. We show 
that the proposed improvements on the model are essential for studying isospin effects in nuclear 
multifragmentation. In particular, these calculations show that accurate comparisons to experi- 
mental data require that the nuclear masses, free energies and secondary decay must be handled 
with higher precision than many current models accord. 
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I. INTRODUCTION 



Experiments have demonstrated that appropriately excited nuclear systems will undergo 
a multifragment disintegration leading to a final state composed of a mixture of fragments 
of charge 3 < Z < 30 and light particles with Z < 2 ^ . Fragments are produced with large 
multiplicities in central heavy ion collisions at incident energies of Ebeam/^ < 100 MeV 
R ^5 in larger impact parameter heavy ion collisions at Ebeam/^ > 200 MeV 0, Q 
and in central light ion induced reactions at Eteam ^5 GeV [7]. Analyses of two fragment 



correlations indicate breakup timescales for these systems that are consistent with bulk 



disintegration 



14 



satisfying an important premise of equilibrium models Jl 



15| that relate multifragmentation to the nuclear liquid-gas phase transition 




Successful comparisons of such models have been made to the measured fragment mul 



tiplicities, charge and energy distributions 



BBQQ. 



Such success, even for reactions 



where a significant collective energy of expansion is observed 0|, implies that these reac 



tions populate a significant fraction of the availab 
such as excited state and isotopic thermometers 



fiphase space. Experimental observables 



20, 



2l|, and the isospin dependence 



22 



of multifragmentation, suggest a degree o: 



energies or smaller systems or both |23 . 



thermalization less complete for higher incident 



25l^a. Such tests, however, have been rendered 



less conclusive by the inability of many current equilibrium models to accurately describe 
the later stages of the breakup where nuclear structure details determine the spectrum of 
excited states and their decay branching ratios. 

Over the years, different versions of the statistical multifragmentation models have been 



developed jl4, 



g8|,|29|,|30|. ased our model upon many of the theoretical 

foundations described in ref. j27| and included the algorithm on partitioning a finite system 
with two components as described in ref. [2^. We call this earlier model SMM85. In 
the improved Statistical Multifragmentation Model called ISMM, we depart from the latter 
approach that the Helmholtz free energies are calculated by carefully including the measured 
states of the fragments, with empirical binding energies and spins [sil 
expressions for these free energies that approach the free energies of refs. 



27 



35,l33|. We obtain 



at excitation 



energies typical of excited multifragmenting systems. The main differences between the 
properties of the hot systems we calculate and those calculated in SMM85 can be attributed 
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to the more accurate expression for the binding energies that we employ; the structure of the 
low-lying states of the fragments plays little role in properties of the hot system. However, 
these structure effects become critical when the fragments cool later by secondary decay. 

Comparisons between results from ISMM and SMM85, reveal large differences between 
the predicted observables, calling many of the previous conclusions into question. In par- 
ticular, we have found that SMM85 calculations tended to overpredict the yields of heavy 
fragments, and consequently, to underestimate those of the lighter ones. More importantly, 
we find that isotopic yields and observables like the isotopic temperatures require careful 
attention to the structure of the excited fragments. If such structural effects are included 
many experimental trends of these observables can be reproduced, and when they are not, 
the experimental and theoretical trends are very different from each other. 

In the following, we recapitulate briefly the formalism of SMM85 and describe in detail 
how we incorporate the improved structure information in the calculation of the properties 
of the hot system at freezeout. This is followed by a description of the secondary decay of 
the hot fragments. Then, we turn to the comparisons of ISMM to predictions of SMM85 
calculations that take less care with these nuclear structure effects. We then compare the 
present improved model to the available experimental data. Finally, we summarize our work 
and provide an outlook towards future comparisons of data to equilibrium models. 



II. THE STATISTICAL MULTIFRAGMENTATION MODEL 



During the later stages of an energetic nuclear collision, the excited system may expand 
to subnuclear density. This expansion may reflect the relaxation of a compressed system 



or the thermal expansion 
3f 



37, or in 



formed in central collisions between comparable mass nuclei 
of a highly excited system formed in a peripheral heavy ion collision 
a collision between a light projectile and a heavy nucleus j^]. For appropriate conditions, 
the excited system disassembles over a time scale of 50-150 fm/c 0, M, till; ll^ into a 
mixture of nucleons, light particles with A< 4 and heavier fragments. Equilibrium models 
[13! such as the statistical multifragmentation models assume that phase space is 

sufficiently well occupied so that the system can be approximated by an equilibrated breakup 
condition characterized by the thermal excitation energy E*, the density p, the mass Aq and 
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the atomic number Zq. Then a second "freeze-out" approximation is invoked, which assumes 
that the system disassembles sufficiently rapidly that further interactions between the various 
particles in the equilibrated breakup can be neglected and that subsequent secondary decay 
of the excited fragments can be calculated as if these fragments are isolated. 

The values of the three conserved quantities E*,Ao, and Zq strongly reflect the dynamics 
of the excitation process and as this dynamics lies outside SMM, they become constraints 
that are introduced as input parameters to the model. The SMM then performs the two es- 
sential tasks required of equilibrium statistical multifragmentation models: (1) the sampling 
of the equilibrium multiparticle phase space, and (2) the secondary decay of excited frag- 
ments. The ffist step in sampling the multiparticle phase space within the SMM is to select 
a fragmentation mode 'm' characterized by a set of particles {A^A,z}m; which are present in 
the equilibrium stage. For each fragmentation mode, mass and charge conservation provides 
that: 

Ao=J2Na,zA and = iV^.z Z (1) 

{A,Z} {A,Z} 

where iV4 ^ is the multiplicity of a fragment, whose mass and atomic numbers are, respec- 
tively, A and Z. The total multiplicity Mm of the fragmentation mode is related to A^a.z 
by: 

Mm= J2 ^^'^ ■ (2) 
{A,Z} 

The selection of the fragmentation modes and the sets of particles {A^_A,z}m for each mode, 
is performed by an algorithm, described in ref. that ensures that all probable choices 
{NA,z}m are sampled, but the frequencies of sampling for the various modes do not re- 
flect their relative contributions to the multifragmentation phase space. This requires the 
introduction of weights Um discussed below. 

The phase space of states consistent with a decay mode {A^A.zjm reflects the number of 
states and consequently the entropy consistent with that mode. The major contributions 
to the total entropy are the entropies corresponding to the internal motion, i.e. internal 
excitation of the fragments. These entropies are calculated within the SMM by introducing 
a temperature Tm for the decay mode. The ensemble average of the expression for energy 
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conservation can then be used to determine the appropriate value of Tm as follows: 

El%, +E* = Ccj^ (^) + J2 N,^,E,M • (3) 

^^^^^ 

On the left hand side of the equation, the total energy is decomposed into the total ground 
state and total excitation energies of the source. The ground state energy, E^-^^^, represents 
the ground state energy of the source calculated as a single spherical nucleus. The first term 
on the right is the Coulomb energy of a homogeneous sphere of matter containing the total 
charge Zq and mass Aq, which is evaluated at a density p = po {Vo/Vmo) where po is the 
saturation density and Vq and Vmo are the volumes occupied by the system at saturation and 
at the breakup densities, respectively. The remaining terms on the right hand side are energy 
contributions, i.e. the kinetic, ground state, extra Coulomb, and excitation energies of the 
individual fragments that are specified below. For the Coulomb energy, this decomposition 
is enabled by invoking a modified Wigner Seitz approximation |4y|, whose accuracy for 



the multifragmentation process has been explored in refs. [221, |2^. The result of applying 
Eq. Q is to obtain values for which conserve energy for the ensemble averaged mean for 
each decay mode and consequently fluctuate from one decay mode to another reflecting the 
corresponding variations in the Coulomb, kinetic and ground state energies of the collection 
of fragments that characterize each decay mode. 

The weight Um for each decay mode is calculated by evaluating the corresponding number 
of states for the mode 

ujm = exp{Sm), (4) 

where the total entropy Sm of the mode is obtained by summing the contributions from each 
particle 



Sm — Sa,z, (5) 



where both E^ z and Sa,z ai'e obtained from the Helmholtz free energies F^^z via the usual 
thermodynamical relations: 

OF 
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and 



F = E-TS , (7) 
which apply to both the contributions from individual fragments and to their overall sums 



of Em, Sm and 



The contributions Fyi z{Tm), associated with each fragment in the partition may be de- 
composed into four terms: 



FaAT) = E^z - Cc ^ 



+ Fl,{T) + FX^T) , (8) 



where -E^'f' is the ground state energy of the fragment. The kinetic term corresponds to: 

'ttinATY''^' 



F^JT) = -Tin 



+ T\n{NA,z^)/NA,z. (9) 

In this expression, V/o = Kno — is the free volume, mj^ represents the nucleon mass, and 
gA,z is the spin degeneracy factor. Empirical ground state spin degeneracy factors are used 
for y4 < 5 because these nuclei have no low lying excited states. For simplicity, we take 
qaz = 1 for heavier nuclei because the influence of non-zero spins on F^z{T) is small and 
can be compensated by small changes in the level density expression for the fragment. The 
Coulomb term, —Cc-^ \ ® represents the corrections in the Wigner-Seitz 

approximation for the individual particles. The excitation of the intrinsic degrees of freedom 
is taken into account by -F|^(T), and is zero for light particles with no excited states. 

To calculate the properties of the multifragment emission from the excited source, one 
should sum the contributions of all the partitions consistent with energy, mass and charge 
conservation. Such a procedure, however, would be extremely time consuming owing to the 
huge number of possible modes. Therefore, the present approach samples the more probable 
modes via a Monte Carlo calculation. This is discussed in detail in ref. {2^; we note in 
passing that the Monte Carlo procedure introduces a bias since not all the mass and charge 
partitions enter with the same weight. Therefore Um must be modified to correct for this 
bias ^]. 



Taking these modifications into account, the average value of a physical observable O is 
calculated by taking a weighted average, 

(O) = ^"^^"^^"^ . (10) 

This average applies both to observables calculated from the primary distributions and 
from the secondary distributions. Because the weights are not unity, the calculation of the 
statistical uncertainties associated with the Monte Carlo procedure requires care. They 
can be easily obtained, however, by repeating the Monte Carlo procedure with a different 
initialization of the random number generator and calculating the variance of the fluctuations 
in the predicted observables. 



A. Ground state energies 

Since the predicted primary yields of excited fragments are exponentially related to their 



binding energies l4| , it is natural to assume that accurate values for the ground state masses 
for the observed fragments are needed. In addition, the isospin dependence of the masses 
and consequent yields of heavier nuclei away from the valley of stability can influence the 
predicted yields of measured light nuclei closer to the valley of stability because all yields 
must be consistent with the constraints imposed by mass and charge conservation. To 
provide more accurate predictions of isotopic distributions, it is relevant to replace the 
somewhat inaccurate Liquid Drop Mass (LDM) parameterization used by many 

current SMM codes QQ- 

To address this problem, we use the recommended binding energies values from Audi 
and Wasptra j^^l when available. The sampling of the most probable partitions discards 
extremely exotic fragments, which would contribute with a vanishing statistical weight. 
Nonetheless, applications of the SMM to realistic multifragmentation scenarios require the 
code to predict the binding energies for many nuclei that have not been measured. Therefore, 
we use a more accurate description of unknown masses given in ref. j^]: 



AV3 



+ 5a,zA~^'^ + , (11) 
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where 



k 



(12) 



A 

and i = V,S stand for volume and surface, respectively. The coefficient 5a,z corresponds to 
the usual pairing term: 

f 

+Spairing N and Z cveu 
Sa,z = < A odd (13) 

^ -^pairing N and Z odd 

The parameters corresponding to the best fit of the empirical masses in ref. are 
ay=15.6658 MeV, ag = 18.9952 MeV, fc=l. 77441, Cc = 0.720531 MeV, 5pairing = 10.8567 
MeV and Cd = 1.74859 MeV. To illustrate the improvement in the model, the upper panel 
(a) of Figure ^ shows the difference between the calculated binding energies from the pa- 



rameterization of the LDM of ref. 28| used in most current SMM codes and the empirical 
values. The lower panel (b) shows the corresponding comparison between the calculated 
binding energies using Eq. ()11|) with the improved parameters (ILDM) and the empirical 
values. One should note that the total binding energies are plotted, rather than the binding 
energy per nucleon. This improved agreement suggests that the predictions for unmeasured 
masses will also be improved. 

Despite the improvement in the overall mass predictions, there can be discontinuity be- 
tween the extrapolated (dashed line) and empirical values (points) as illustrated in Fig. |21 
To improve the matching between the binding energies of the known masses and the ones 
predicted by our mass formula, we compute average shifts of the ILDM formula from the 
empirical values and use these shifts to correct the values in Eq. (fTTj) . For an isotone that 
has a lower charge than its isotonic partners in the compilation of ref. js^ we use the three 
lightest isotones with the same value of N in the compilation to compute the shift. Simi- 
larly for an isotone that has a higher charge than its isotonic partners in the compilation 
of ref. we use the three heaviest isotones in the compilation to compute the shift. This 
shift is then subtracted from the prediction of the ILDM formula: 

jjextrap _ r,ILDM \ /-,a\ 
^A,Z - ^A,Z - , (14) 
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where 



A«=5i:Kt?f'-<"t.). (15) 

i=l 



^Aj^Zi is corresponding value from the compilation of ref. [3^. Two shift values are 
therefore computed for each value of N. The final binding energy values used in the ISMM 
calculations are illustrated for four cases by the solid lines in Fig. |21 where it shows that the 
discontinuity between the empirical (star) and extrapolated (dashed line) values is removed. 



B. Fragment internal free energies 

In this work, we have modified SMM85 so as to allow accurate predictions of isotopic 
properties, but have limited the extent of these modifications in an effort to retain many of 
the predictions of the original theory. In particular, we have retained the high temperature 
properties of the fragment free energies, ^, which are parameterized here and in the 
SMM85 as: 



rp2 _ 2^2\ 5/4 



rp2 

- A- , (16) 

Co 



where jSo = 18.0 MeV, eo = 16.0 MeV, and Tc = 18.0 MeV. This expression holds only 
for temperatures smaller than critical temperature, Tq- At low temperatures, T « Tc, 
this expression depends quadratically on T as expected for a Fermi liquid. At the critical 
temperature where the surface tension vanishes, the surface energy contribution to the total 
free energy Fa,z{T) falls to zero when the surface energy contribution in Eq. ()16|1 is combined 
with the corresponding ground state energy term in Eq. (jHl) • As we do not calculate decays 
at T > 10 MeV, we do not concern ourselves here with the form for F|^(T) at T > Tq- 
For 3 MeV ^ T < 10 MeV, where multifragmentation is important, however, this form for 
FXz{T) in Eq. ()16|) is not unique, and other expressions with different thermal properties 
should be explored. In the following we introduce empirical modifications to this free energy 
expression by taking into account the nuclear structural information of known excited states. 

First we turn our attention to the fact that most fragments at T > 2 MeV are parti- 
cle unstable and will sequentially decay after freezeout. This decay is sensitive to nuclear 
structure properties of the excited fragments such as their nuclear levels, binding energies, 
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spins, parities and decay branching ratios. The first three of these quantities also influ- 
ence the free energies; this can be calculated via the fragment internal partition functions. 
Self-consistency in the freeze-out approximation dictates that the states from which these 
fragments decay after freezeout should be consistent with the Helmholtz free energies used 
in calculating the primary yields of the hot fragments at freeze-out. 

In order to discuss this self-consistency requirement, we must consider the density of 
states pstates{E) and its mathematical relationship with the Helmholtz free excitation energy 



where the integral is over the excitation energy E of the nucleus. Here we have, for simplicity, 
neglected the complications of a degenerate ground state, which contributes negligibly to the 
free energy at high excitation energy. In the original papers on the SMM, the level densities 
corresponding to the SMM were not stipulated. We now consider what is required of the 
density of states to achieve the high temperature behavior for F^^(T) given by Eq. (jlfjj] . 
Then we will address the general issue of making the level densities consistent with empirical 
information and how that impacts the free energies. Finally, we will discuss specific details 
of the incorporation of the empirical information into the level density expressions. 

1. High temperature behavior 

First we investigate what forms of level densities may be consistent with the free energies 
in Eq. (fT^ . We note that the functional dependence of F|^(T) used in Eq. lfTBjl makes its 
analytical inversion difficult at high temperatures. Instead, it is easier to find a smooth 
real functional form for pstates{E) that reproduces the numerical values for FX z(T) at high 
temperatures than it would be to perform an inverse Laplace transformation of F^ziT) in 
the complex plane. We note that if one inverts a Taylor expansion of F|^(T) up to second 
order in T by the saddle point approximation, one obtains the Fermi gas expression: 



where asuM is the absolute value of the coefficient of the second order term of the free 



F*{T): 




(17) 




1/4 



(18) 
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energy expansion in T: 



asMM = — + 7;Po^^- (19) 

However, this expression is unsatisfactory at high temperatures, as is illustrated in Fig. El 
when the free energies obtained from Eq. ^1 (dashed lines) are compared with SMM85 free 
energies in Eq. UHl (solid lines). Instead, we take Eq. (THlas a starting point and obtain a 
useful analytic expression by multiplying PFG,states{E) by an ad hoc energy dependent term 
to obtain free energy values in numerical agreement with Eq. ()16|) : 

PSMM,states{E) = PFG,states{E) e-^SM«(»SMMi?)^/^ ^ (20) 

where bsuM is given by: 



hsMM = 0.07A-^ (21) 

r = 1.82 (1 + . (22) 

V 4500y ^ ^ 

The free energies obtained via Eqs. (|17p and (|2Uj) are displayed in Fig. |21 as symbols for 
two different mass regions. This simple parameterization is fairly accurate at temperatures 
T < 10 MeV, in the range of interest. 

2. Empirical Level densities at low excitation energies for Z<15 

Several factors motivate the efforts to develop an accurate treatment for the level densities 
at low excitation energy for Z < 15. The first factor is that most multifragmentation data 
are available for light fragments in this mass range. The second is that empirical nuclear 
structure information is also available for these nuclei. A comparable treatment of the level 
density for the heavier fragments would be interesting, but the needed structure information 
is frequently incomplete or entirely missing. Fortunately, if we focus on the yields for A < 8, 
the contributions from the secondary decay of the heavy nuclei with Z > 15 are of the order 
of 10%. Thus the errors introduced by the neglect of this structure information for the heavy 
nuclei does not strongly influence the results of the final yields and one can proceed towards 
reasonable predictions at the present time. 
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At lower excitation energies, it is customary to discuss the density of levels pieveis rather 
than the density of states because this definition is more useful experimentally when the 
spins of specific levels are not accurately known. Mathematically, the density of states is 
related to the densities of levels for individual spin values pieveis{E, J) by: 

Pstates (^) = 5^ (2 J + 1) pievels (E, J) (23) 
J 

While the spacings between energy levels in a given nucleus generally decrease smoothly 
with excitation energy, as a practical matter one often decomposes the empirical level den- 
sity Pemp,ieveis{E, J) iuto two cxpressious that apply in two different approximate excitation 
energy domains: (1) one (labelled as PD,ieveis{E, J)) containing discrete well separated states 
at low excitation energies and (2) another (labelled as pc,ieveis{E, J)) containing a continuum 
of overlapping states at higher excitation energies. For Z<15, empirical level information 
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33l | is applied as much as possible to the low-lying discrete level density, wherever the 



experimental level scheme seems complete, 

PD,leyels {E,J) = Y,m-E), (24) 

i 

where the summation runs over the excitation energies Ei corresponding to states of spin 
J . Examples of empirical levels for ^°Ne and '^^P are shown as bars in Fig. |3] For higher 
excitation energies, a good approximation to the continuum level density has been obtained 
by ref. 4^ by combining Fermi liquid theory, a simple spin dependence and experimental 
knowledge. The relevant expressions, shown as dashed lines in Fig. HI are 3| , 

Pc,ieveis{E, J) = pc{E)f{J, a) for E>E, (25) 

where 

_ exp[2v/a(E - Ep)] 
'^^""^ - l2V2ay\E - Epfl^a ^'^^ 
, . (2J+l)exp[-(J+l/2)V2a^] 

(7^ = 0.0888^/a{E - Eo)A'^/^, (28) 



and the level density parameter a = A/8. E*, J, A and Z are the excitation energy, spin, 
mass and charge numbers of the fragment. Eq is determined by matching the total high-lying 
level density to the total low-lying level density as follows, 

i>Ec r pEc p 

j dE dJpc,levels{E,J) = / dE dJpD,levels{E,J), (29) 

Jeo J Jo J 
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where Ec is the energy at which the switch from discrete to continuum level density expres- 
sions is made. 

The comparison in Eq. (j29|) is between the total level densities summed over spin. This 
is done primarily to reduce the sensitivity in the matching to uncertainties in the spin 
assignments for some of the discrete states. By adjusting the parameter Eq, the total level 
density for continuum states was connected smoothly to the total level density for low-lying 
states a.t E < Ec and Z < 12. The connection point Ec to high- lying states, for Z < 12, 
was chosen to be the maximum excitation energy up to which information concerning the 
number and locations of discrete states appears to be complete so that the empirical level 
density fEql24|) was solely applied for low-lying states. 

For the case of Z > 12, low-lying states are not well identified experimentally and a con- 
tinuum approximation to the discrete level density 3| was used by modifying the empirical 
interpolation formula of Ref. 4^ to include a spin dependence: 

PD,ieveis{E,J) = ^ exp[{E - Ei) /Ti] 
1 



^ (2J + l)exp[-(J + l/2)V2ag] 

E.(2J. + l)exp[-( J, + 1/2)72^0^,]' ^ > 



for E < Ec, where the spin cutoff parameter al = 0.0888^/a{Ec - Eo)A^/^. For Z > 12, the 
values of Ec = Ec{A, Z) were taken from Ref. J45|] as well as parameters Ti = Ti{A, Z) and 
El = Ei{A, Z), and in this case, the approximate level density fEq]30|) was used in place of 
an empirical level density for the low-lying states. 



3. Matching low and high excitation energy behavior 

Now, we turn to the requirement of self consistency between the expression for ^^^(T) 
and the level density relevant to secondary decay. In general, secondary decay becomes more 
sensitive to nuclear structure quantities such as the excitation energies, spins, etc. as the 
systems decay towards the ground state. At low excitation energies, one is more accurate 
using empirical level densities in place of the expression in Eq. ( I16|) . which does not even 
depend on Z. As the excitation energy is increased, however, the continuum level density 
becomes very large, little sensitivity to nuclear structure details remains and a simpler 
expression like Eq. (fTHI) may suffice. 
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In the following, we take psMM,states{E) to be the state density at high energies and 
match it to the continuum part of the empirical state densities at low excitation energies. 
This procedure uses the empirical information for excitation energies E* < E^, a linear 
interpolation for < E* < + AE, and psMM,states{E) at higher values of the excitation 
energy. The net result is a set of level density and state density expressions that span the 
range of excitation energies relevant to multifragmentation phenomena. For E* < Ec, one 
uses the expression for the discrete, low-lying state density, 

PiSMMiE,J)=pniE*,J). (31) 

For Ec < E* < Ec + AE, the new level density is an interpolation involving the continuum 
expression relevant at low excitation energies between pc, states and Psmm, states, 

E* -Er. 



Pismm{E*,J) = Pc{E*,J){l 



AE ' 

+Psmm{E\J) ^*~^\ (32) 

where AE = 2.5A MeV provides a smooth transition from pc to Psmm- The SMM level 
density (shown as dotted lines in Fig. can be incorporated with a similar spin dependence 
as in Eq. ESI 

PsMMiE*,J)=psMMiE*)f{J,a). (33) 

For E* > Ec + AE, the new density simply becomes the same as the SMM level density 

Psmm, 

PiSMM — Psmm {E*,J). (34) 

In Fig. m the empirically modified level density described in Eqs. (I3HI34|1 is plotted as solid 
lines for ^°Ne and ^^P. 

The level density pc in Eq. ()25|) can be used as a proper extension to the low-lying level 
density p_D in Eqs. and (IHUj) and a bridge for matching to the SMM level density at 
continuum. Such a matching procedure provides a state density that is empirically based at 
low excitation energies but becomes progressively more uncertain as the excitation energy 
is increased above E* ^ Ec- This uncertainty in the thermal properties of nuclei at such 
high excitation energies is not a question of finding an appropriate interpolation, but is, in 
fact, a fundamental issue that must be resolved by comparisons to experimental data. For 
example, other expressions can be proposed for the level density at E* > Ec and comparisons 
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of experimental data to SMM predictions of sensitive multifragment observables can be used 
to constrain the level densities at high excitation energies. 

Free energies F^^{T), which reflect contributions from the discrete excited states are 
obtained by inserting this parameterization for Pstates{E) into Eq. (jl7|) . and performing 
a numerical integration. To facilitate the insertion of these free energies into the SMM 
algorithm, we parameterize ^^^(T) by: 

Ea,z(T) = Fgj^.j^^{T) ( 1 - — -— — -— — — r ) , (35) 

where Fgj^j^j{T) stands for the SMM internal free energy of Eq. (fT?)|l which is adopted in 
various SMM models. The parameters Tadj and AT are adjusted to reproduce the numerical 
calculation of F*{T) provided by Eqs. (Hlj) and (PTTllHij) for T < 10 MeV. In these fits, a 
value for AT = 1.0 MeV is used for most nuclei (The exceptions are mainly very light 
nuclei.), while Tadj is varied freely. The accuracy of the fit is illustrated in Fig. which 
compares the exact values of F*{T) (symbols) to the approximation given by Eq. (jHSj) (solid 
line), for a ^''Ne nucleus. The dashed line in this figure represents the free energy used in 
SMM85 calculations in which the experimental discrete levels are neglected. The matching 
procedure allows the discrete excited states to dominate the low temperature behavior, while 
the high temperature behavior remains similar to that of the SMM85, consistent with the 
goals stated above. 

Because the empirical level densities vary from nucleus to nucleus, the parameters Tadj 
and AT must be fitted for each nucleus used to obtain T^^(T). Fits of the same quality 
as that for ^''Ne are achieved for all the light nuclei with Z < 15. These fitted values of 
Tadj are shown as symbols in Fig. IHl We do not perform such fits for Z > 15, because the 
level density information there is less complete. We nevertheless extrapolate the main trend 
of the parameters to heavy nuclei, for which detailed experimental information on discrete 
excited states is not available, in order to avoid spurious discontinuities in the equilibrium 
primary yields. As mentioned above, there seems to be a very weak dependence on AT and, 
therefore, we assume AT = 1.0 MeV for Z > 15. In spite of the uncertainty in extrapolating 
Tadj , the dashed line in Fig. IHl shows that 

Tadj = 22.0^-°-^ MeV {Z > 15) (36) 
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describes the trend (dashed hne) for the lower masses and we adopted it for the higher 
masses as weU. 

III. SECONDARY DECAY 

With few exceptions, the stable yields after secondary decay are the quantities that are 
usually measured experimentally. An accurate secondary decay procedure is indispensable 
to calculate the contributions from secondary decay and deduce the information of the 
primary hot system from experimental data. The sequential decay procedure consists of 
two parts. One is to decay particles with Z<15 through a large empirical (MSU-DECAY) 
table including all the states of nuclei with known information such as binding energy, s pin, 
isospin, parity and decay branching ratios. The other part is to use the Gemini code j47| 
for particles outside the empirical table (usually Z>15). 

A. Decay table 

The implementation of Eqs. (j31H34j) involves the construction of a 'table' of quantities 
such as the excitation energies, spins, isospins, and parities of levels of nuclei with Z < 15. 
For excitation energies E < and Z < 15, each of the entries in the table corresponds 
to one of the tabulated empirical levels. When the information on the level is complete, it 
is used. For known levels with incomplete spectroscopic information, values for the spin, 
isospin, and parity were chosen randomly as follows: spins of 0-4 (1/2-9/2) were assumed 
with equal probability for even-A (odd-A) nuclei, parities were assumed to be odd or even 
with equal probability, and isospins were assumed to be the same as the isospin of the 
ground state. This simple assumption turns out to be sufficient since most of spectroscopic 
information is known for these low-lying states. 

For excitation energies where little or no structure information exists, the level density 
was assumed to be given by the level density algorithm discussed in the previous section 
and groups of levels were binned together in discrete excitation energy intervals of 1 MeV 
for E* < 15 MeV, 2 MeV for 15 < < 30 MeV, and 3 MeV for E* > 30 MeV in order to 
reduce the computer memory requirements. The results of the calculations do not appear to 
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be sensitive to these binning widths. A cutoff energy of E*^^^jrj/A = 5 MeV was introduced 
corresponding to a mean hfetime of the continuum states at the cutoff energy about 125 
fm/c. For simphcity, parities of these states were chosen to be positive and negative with 
equal probabihty and isospins were taken to be equal to the isospin of the ground state of 
the same nucleus. 



B. Sequential decay algorithm 

Before sequential decay starts, hot fragments from primary breakup need to be populated 
over the sampled levels in the prepared table according to the temperature. For the ith level 
of a given nucleus (A,Z) with its energy E* and spin Jj, the initial population is, 

v^ (2J. + l)exp(-E;/r)p(E;,J,) 
' ~ ' ^ E.(2 J. + 1) eM-E:/T)p{E:, J.) ^^'^ 

where Yq is the primary yield of nucleus (A,Z) and T is the temperature associated with the 
intrinsic excitation of the fragmenting system at breakup. 

Finally all the fragments will decay sequentially through various excited states of lighter 
nuclei down to the ground states of the daughter decay products. Eight decay branches 
of n, 2n, p, 2p, d, t, ^He and alpha were considered for the particle unstable decays of 
nuclei with Z< 15. The decays of particle stable excited states via gamma rays were also 
taken into account for the sequential decay process and for the calculation of the final 
ground state yields. If known, tabulated branching ratios were used to describe the decay 
of particle unstable states. Where such information was not available, the branching ratios 

n 

were calculated from the Hauser-Feshbach formula |48l |. 
where 

\Jd + Je\ \Jp + J\ , , _,w 

X E E '^''T (39) 

J=\Jd~Je\ l=\Jp — J\ 

for a given decay channel d (or a given state of the daughter fragment). Jp, J^, and Je 
are the spins of the parent, daughter and emitted nuclei; J and / are the spin and orbital 
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angular momentum of the decay channel; Ti{E) is the transmission coefficient for the Ith 
partial wave. The factor [l + 7rp7rd7re(— 1)']/2 enforces parity conservation and depends on the 
parities vr = ±1 of the parent, daughter and emitted nuclei. The Clebsch-Gordon coefficient 
involving Ip, 1^, and Je, the isospins of the parent, daughter and emitted nuclei, likewise 
allows one to take isospin conservation into account. 

For decays from empirical discrete states and / < 20, the transmission coefficients were 
interpolated from a set of calculated optical model transmission coefficients; otherwise a 



parameterization described in Ref. 



was applied. 



IV. MODEL PREDICTIONS AND COMPARISONS 
A. Caloric Curve 

Before presenting predictions for isotope distributions and other observables for which 
the present theoretical developments were undertaken, we examine predictions of the present 
improved model for the caloric curve and the primary fragment multiplicities, both of which 
displayed features in SMM85 and other SMM calculations ^ that are characteristic of 
low density phase transition. For example, SMM85 calculations predict an enhanced heat 
capacity for multifragmenting systems reflecting the latent heat for transforming nuclear 
fragments (Fermi liquid) into nucleonic gas. FigEI shows the caloric curve, i.e. the depen- 
dence of the mean fragmentation temperature (T^) on excitation energy, for a system with 



Ao=168 and Zo=75. In 
by the original SMM85 



3oth panels, the dotted lines indicate the relationships predicted 
271 |28|, the solid lines indicate the corresponding predictions of 
the ISMM with all the modifications discussed in this paper and the dashed lines indicate 
the results provided by an SMM85 calculation that uses the new binding energies of Eqs. 
()lHll5j) and the old parameterization of ref. j^] for the Helmholtz free energies. These 
latter calculations allow one to assess the impact of the changes in the binding energies and 
free energies independently. 

The two panels provide the caloric curves corresponding to two different constraints on 



the density. In the lower panel, a multiphcity-dependent breakup density is assumed, 
corresponding to a fixed interfragment spacing at breakup; this leads to a pronounced plateau 
in the caloric curve for all three calculations. By taking into account the kinetic motion and 
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the Coulomb interaction, we have estimated the pressure using the relationship 



where P is the pressure, M is the total multiplicity, Vf is the free breakup volume and V 
is the total volume. Limiting the pressure estimates to temperatures for which the multi- 
plicity exceeds 10 and the pressure can be more reasonably defined, we show the pressure 
corresponding to these multiplicity-dependent breakup densities in the lower panel of Fig. 
IHl The corresponding primary fragment multiplicities are shown in the lower panel of Fig. 



nultipl 

• 3, 



El Consistent with the conclusions of ref. |49|, we find the requirement of approximately 
constant interfragment spacing corresponds to breakup pressures that exhibit only a small 
fractional increase with temperature. In the upper panels, we show the corresponding caloric 
curves (FigEj), pressures (Fig. IH]), and multiplicities (Fig. ^ calculated at fixed breakup 
density p/po = 1/6. These show a steeper dependence of the caloric curves on excitation 
energy and the small maximum displayed in the lower panel of Fig|7| at excitation energies 
of about 3 MeV disappears. The corresponding pressures at constant density, shown in the 
upper panel of Fig. |H1 increase monotonically with excitation energy. However, they are 
lower than those calculated assuming a multiplicity dependent breakup density, because the 
density for the constant volume calculations is lower. 

These figures reveal that the trends of the thermal dynamical properties of these three 
models to be similar. In general, the temperatures in the plateau region atE*/A = 3—8 MeV 
in the lower panel of Fig. 7 are larger for the ISMM calculations using the improved free 
excitation energies. This is consistent with the fact that the level densities and, consequently, 
the entropies of the fragments are lower in the improved model, which generally raises the 
temperature corresponding to a given excitation energy. Specifically in the plateau region, 
reducing the entropies of the fragments raises the latent heat for the transformation from 
excited fragments to nucleon gas and raises the temperature at which the transition occurs. 
The influence of the improved binding energies on the caloric curve is less obvious, but this 
change seems to be largely responsible for the differences between the SMM85 and ISMM 
at E*/A>6 MeV. 

Discussions of the nuclear caloric curve usually focus on the excitation energy dependence 
of the temperature and ignore the density dependence. To illustrate that the phase diagram 
is two dimensional and a density dependence does exist, we contrast in Fig. [TUlthe density 
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dependence (right panel) of the temperature at a fixed excitation energy of E*/A = 6MeV 
(open squares) to the excitation energy dependence (left panel) of the temperature at a 
fixed density of p/po = 1/6 (solid circles). Both the excitation energy and the density 
dependences of the caloric curve are clearly important. It is therefore relevant to find and 
measure observables that constrain significantly the freezeout density. 

B. Charge and Mass Distributions 

Calculations of the mass distribution (left panel) and charge distribution (right panel) for 
excited primary fragments are shown in Fig. ^2 for a system with Aq = 186 and Zq = 75 at 
E*/A = 6MeV. This system has the same charge to mass ratio as the symmetric ^^^Sn+^^^Sn 
system, but is chosen to be 3/4 of the total mass in order to approximately address the mass 
loss to preequilibrium emission. The dotted lines denote the predictions using SMM85 and 
the solid lines denote the predictions using ISMM. The primary distributions from ISMM 
fiuctuate about the smooth distributions of SMM85 for Z < 20 and A < 60 and then fall 
below SMM85 at higher mass and charge. The fiuctuations are related to the infiuence 
of shell and pairing effects on the ground state masses, which have no significant impact 
on the final yields after secondary decay as discussed below. The trend of reduced yields 
at higher masses and charges is related to the tendency shown in Fig. ^ for the binding 
energies in the SMM85 to consistently exceed the empirical values at Z > 20 and A > 60. 
Because conservation of mass and charge dictates that an increase in the yields of heavier 
fragments must be compensated by a decrease in the yields of the lighter ones, one should 
see a comparable under-prediction of the primary yields of the lighter fragments by SMM85. 

Fig. ^1 shows the corresponding final mass (left panel) and charge (right panel) distri- 
butions after secondary decay. The solid lines denote the predictions using ISMM. Experi- 
mental fragment yields from the central ^^^S'n +^^^ Sn collisions are plotted as solid points 
Isoj ]. To investigate the infiuence of the fiuctuations in the primary distributions due to 
shell and pairing effects on the ground state masses, we have decayed the primary fragments 
from the SMM85 via the same empirical secondary decay procedure discussed in Sect. IIIII 
The final mass and charge distributions of the SMM85 are shown as the dashed lines in 
Fig. El Minimal discrepancies are seen in low mass and charge regime indicating that the 
secondary decay mechanism washes out the fiuctuations in the primary distributions due to 
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the influence of sliell and pairing effects on tlie ground state masses. Meanwliile, significant 
differences on lieavy fragments remain. In order to see tfie differences between tfie two cal- 
culations, the low A and Z regions are expanded in Figure E| Here again, experimental 
fragment yields from the central ^^^5*77, +^^^ Sn collisions are plotted as solid points [50.] . 
The agreement is very good even though no special attempt has been made to optimize the 
parameters of the calculations to achieve the best representation of the data. 



C. Isotopic distributions 

In Fig. El the primary isotopic distributions for four elements emitted are shown for a 
system with Aq = 186 and Zq = 75 aX E* /A = 6 MeV. The solid lines show predictions for 
the present improved model and the dashed lines show predictions of the SMM85 code of refs. 



l28j . The two calculations produce primary isotopic distributions that are considerably 





broader and more neutron rich than corresponding final distributions after secondar y d ecay 
shown in Fig. ^1 For reference, the measured isotopic distributions of refs. 0,13 are 
shown as solid points in Fig. ^1 While the parameters of the code were not optimized 
to reproduce the data, it is interesting to note that the widths of the distributions from 
ISMM calculations and data are similar although the data seem to be more neutron rich 
than the calculations. Studies have shown that the final isotopic distributions calculated 
with an empirical secondary decay procedure such as that employed by the ISMM are much 
broader and more neutron-rich than the corresponding distributions predicted by the more 
schematic statistical models . In order to compare with the available experimental data, 
he isospin observables derived from these isotopic distributions such as isoscaling parameters 
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42| and isotopic temperatures require an accurate secondary decay approach with 
detailed nuclear structure information taken into account. 

Isotope thermometers have been utilized as the primary probes for extracting the caloric 
curve of the nuclear liquid-gas phase transition. Since these observables are constructed from 
the isotopic distributions, they share the sensitivity to structure effects in the secondary 
decay discussed above. In the isotopic thermometer technique, the temperature is extracted 
from a set of four isotopes produced in multifragment breakups as follows , 

"° ln{aR) ^ ^ 
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where 



r(A2,Z2)/r(A2 + i.Z2)' ^ ' 

AB = B{Ai,Z,)-B{Ai + l,Z,) 

-5(^2,^2) +5(^2 + 1,^2), (43) 



and 



(2JZ2,A2 + 1) (2Jzi,Al+l + i: 

a 



^i<^l'^ (44) 

,A2+1 "T i 

Here Y{A,Z) is the yield of a given fragment with mass A and charge Z; B{A,Z) is the 
binding energy of this fragment; and Jz^a is the ground state spin of the nucleus. Although 
this expression is derived within the context of the grand canonical ensemble, it has been 
applied to a wide variety of reactions and regarded as an effective or "apparent" temperature 
that may differ somewhat from the true freezeout temperature T. The relationship between 
Tiso and T can be calculated within an appropriate statistical model for the fragmentation 
process if one exists. In general, one chooses a set of four isotopes with large AB to minimize 
sensitivity to details of the corrections from secondary decay. 

To examine the influence of secondary decay, measured and calculated temperatures are 
extracted from double ratios of Z=2-8 fragments and plotted in Figure ^1 The large AB 
requirement generally limits the apparent temperature observables to three types of ther- 
mometers: a.) TiU^^^He), ^2=2, ^2=3, b.) T,,oC^^^^C), Z^ =6, ^2=11, and c.) Ti^X^^^'^O), 
^2=8, y42=15, where thermometer (a) involves the light particle pair ^■'^He while thermome- 
ters (b) and (c) concern only the intermediate mass fragments (IMF's) of Z=3-8. Table I 
lists the corresponding thermometers plotted in Figure El • The top left panel in Fig. IT6l 
shows the ISMM predictions for these three types of thermometers as a function of Ai. 

Since the denominator in Eq. ()42|1 is fixed by classifying the temperatures into three 
types, the fluctuations are related to Ai. In all cases, the two thermometers involving ^°Be 
and ^^O are much higher than the others due to many low lying states in these nuclei 21 1. 
The extracted temperatures from all the other thermometers are significantly lower than the 
primary temperature of 5 MeV which is shown as the dotted line in the four panels. There 
seems to be a Z dependence in Tiso- ^iso(^^'^^O) is about 0.5 MeV lower than Tisoi^^'^'^C) 
which is only slightly lower ( 0.2 MeV) than Tisoi^'^He). In addition, there is also a trend of 
isotopic temperature values decreasing as a function of Ai. The lower temperatures reflect 
increasing contributions of multi-step secondary decay contributions. As these multi-step 
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contributions originate from the decay from an ensemble of unstable nuclei that are less 
excited than the original ensemble, it has the effect of making the system appear cooler. 

For comparison, we use the corresponding isotope temperatures extracted from the data 
obtained in the central collisions of ^'^^Sn +^^^ Sn reactions at E/A=50 MeV |50] shown 
as solid squares(top right panel), circles (bottom left panel) and stars (bottom right panel) 
for Tisoi^'^^He), Tisd^^^^^C), Tiso{^^^^^0), respectively in Figure CHI The calculated ISMM 
isotopic temperatures (lines) follow the trends of the corresponding experimental values. 
Despite the fact that the parameters in the ISMM calculations have not been optimized, the 
calculated temperatures of TisoC^'^'^C) and Tjso(^^'^^O) (bottom panels) are nearly the same 
as the data within the theoretical uncertainties, which indicates that the IMF's distributions 
can be well reproduced in an appropriate equilibrium model. 

However, the experimental Tiso{^''^He) temperatures (solid squares in top right panel) 
are systematically higher than the corresponding ISMM values (dot dashed line). As these 
thermometers derive their sensitivity to the temperature from the large binding energy 
difference between ^He and ^He, the difficulty in reproducing these quantities may arise 
if there are significant nonequilibrium production mechanisms for light particles such as 



^He m 



52| |. To illustrate this effect, we assumed that 2/3 of the measured ^He yield is 
of a non-thermal origin. This increases the '^He yield by a factor of three and the new 
calculations are shown as the solid line in the top right panel. The resulting apparent 
temperatures are nearly the same as the experimental data. This simple assumption explains 
the discrepancies between Tiso(^'^He) and TisoC'^'^'^C) observed experimentally. However, 
the present calculations also suggest that sequential decays have a much larger effects on 
Ti,o{^^^^^C) and Ti,o(^^^^^0) than previously assumed [24 1. 

To illustrate the importance of using an accurate sequential decay code to decay the 
primary hot fragments before data can be accurately compared. Table I contains the ex- 
perimental measured isotope temperatures in the fourth column. Predicted temperatures 
from the ISMM using the MSU-DECAY code are plotted in the fifth column. As shown in 
Figure 16 and Table I, there is a close correspondence in the fiuctuations of the temperature 
between the ISMM and observed temperatures. However, if one uses the SMM code of Ref. 
[4], which contains a Fermi-break up decay mechanism for excited fragments and utilizes 
schematic structure information to calculate the secondary decays, the fiuctuations in the 
temperature, listed in the last column in Table I, are much larger than those observed in 
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the data. In this respect, one should especially note those involving ^'^Li; ^'^°Be, ^^'^^B, and 
I'^'i^O where the calculated Tapp differ from the data by more than a factor of two. The 
discrepancies in the predicted ratios are significantly larger still, by a factor of AB/Tiso, 
according to Eq. (41). 



V. SUMMARY 

The multifragmentation of excited nuclear systems produces excited fragments that decay 
into the observed ground state nuclei by mechanisms that are strongly influenced by the 
ground and excited state spins and energies of the fragments and by their decay branching 
ratios. Prior equilibrium multifragmentation models employed approximate descriptions for 
these quantities that are insufficiently accurate to describe the new isotopically resolved data 



now becoming available j22,|50|. In this paper, we include this information self-consistently, 
building the improved statistical multifragmentation model (ISMM) upon the foundations 
of ref. j22l,|2S|. The main differences between the properties of the hot systems we calculate 



and those calculated in ref. [27|, |28[ can be attributed to the more accurate expression for 
the binding energies that we employ; the structure of the low-lying states of the fragments 
plays little role in properties of the hot system. These structure effects become critical when 
the fragments cool later by secondary decay. 

Our calculations call many of the previous conclusions of equilibrium multifragmenta- 
tion models into question. In particular, we have found that the SMM85 and other similar 
calculations tended to overpredict the yields of heavy fragments, and, consequently, to un- 
derestimate those of the lighter ones. More importantly, we find that isotopic yields and 
observables like the isotopic temperatures require careful attention to the structure of the ex- 
cited fragments. Thus, prior calculations of these isotopic observables using models that do 
not include such structure information accurately may be unreliable and lead to questionable 
conclusions. 
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FIG. 1: Differences between the total binding energies predicted by mass formulae and those 
recomended in ref. [s^. Upper panel :Plot (a) displays the differences when one uses the LDM 
mass formula used in SMM85 calculations. i28ii. Lower panel: Plot (b) displays the differences 
when one uses the ILDM mass formula presented in this work. 
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TABLE I: List of isotopic thermometers plotted in Figure 16. 

IMF-meters AB a Ti^o (Data) Ti^o (ISMM) Tiso{SMM)[20] 



6,7Li/ll>12c 


11.472 


5.898 


3.740 


3.315 


3.625 


7,8Li/ll,12c 


16.690 


5.361 


3.244 


3.212 


4.419 


8,9Li/ll,12c 


14.658 


3.351 


3.146 


3.065 


1.014 


9,iOBe/ii,i2c 


11.910 


1.028 


5.643 


5.102 


12.561 


ll,12gyll,12Q 


15.352 


3.000 


3.651 


3.154 


3.928 


12,13b/11,12q 


13.844 


5.278 


3.720 


3.031 


1.636 


12,13q/11,12(-. 


13.776 


7.917 


3.418 


3.078 


3.608 


13,14q^11,12q 


10.545 


1.962 


3.288 


2.949 


2.590 


15,16]vj/ll,12(-. 


16.233 


9.669 


2.767 


2.564 


2.716 


16,17q /11,12q 


14.578 


23.069 


2.648 


2.443 


2.555 


17,18o/11,12q 


10.678 


0.637 


6.921 


6.009 


4.514 


6,7Li/15,16o 


8.413 


3.050 


2.273 


2.352 


2.209 


7,8Li/15,16o 


13.631 


2.773 


2.636 


2.565 


3.084 


8,9Li/15,160 


11.599 


1.733 


2.476 


2.368 


0.768 


9,10Bg/15,16o 


8.851 


0.532 


4.143 


3.610 


5.562 


ll,12g^l5,16o 


12.293 


1.551 


2.906 


2.466 


2.701 


12,13b^15,16o 


10.785 


2.729 


2.109 


2.303 


1.184 


12,13q/15,16o 


10.717 


4.094 


2.643 


2.334 


2.402 


13,14q^15,16q 


7.486 


1.014 


2.316 


2.270 


1.588 


15,16j^ ^15,160 


13.174 


5.000 


2.236 


2.043 


1.990 


16,170/15,160 


11.519 


11.930 


2.083 


1.893 


1.814 


17,180/15,160 


7.619 


0.330 


4.863 


4.027 


2.523 


6.7Li/3.4He 


13.328 


2.183 


5.693 


3.632 


4.708 


7'8Li/3.4Hc 


18.546 


1.984 


4.197 


3.431 


5.386 


8'9Li/3'4He 


16.514 


1.240 


4.200 


3.309 


1.169 


9,10Bg/3,4jje 


13.766 


0.380 


9.938 


5.413 


22.410 


ll,12B/3,4jje 


17.208 


1.110 


4.948 


3.390 


4.814 


12,13B/3,4jjg 


15.700 


1.953 


3.599 


3.287 


1.931 


12,13c/3,4jj^ 


15.632 


2.930 


4.731 


3.337 


4.487 


13,14c/3,4jje 


12.401 


0.726 


1^00 


3.276 


3.319 


15,16N/3,4Hg 


18.089 


3.578 


3.519 


2.766 


3.206 
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FIG. 2: Total binding energies for various nuclei. The full lines correspond to the corrected 
LDM formula, whereas the symbols represent the experimental data of ref. [s^. The dashed lines 
correspond to the predictions given by Eq. Ulljl using the optimized parameters. For details see 
the text. 
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FIG. 3: Internal free energies for A = 20 (upper panel) and A = 200 (lower panel). The SMM85 
expression [Eq. (fT^ ] is represented by the full line whereas the dashed lines stand for the results 
obtained with the Taylor expansion [Eq. (|18|) ]. The Free energy calculated through the level density 
given by Eq. (PU)) is depicted by the symbols. 
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FIG. 4: Level densities as a function of excitation energy for ^"Ne and ^^P. Two energy ranges are 
plotted to show the behaviors of level densities at both low (left panels) and high energy (right 
panels) ends. The density of experimentally known levels is shown as bars in the low energy region. 
The dashed lines are the extrapolations of the empirical values according to Eq. 1251 The dotted 
lines are the level density (Eq. I2U|) parametrized from the SMM85. The solid lines are the level 
density adopted in this work (Eqs. I^Tll51]l . 
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T (MeV) 

FIG. 5: Comparison between F*{T) calculated through Eqs. ()17() and H31() - H34() . symbols, and the 
approximation given by Eq. H35|) . full line. To illustrate the influence of quantum effects at low 
temperatures, the dashed line represents the free energy used in SMM85 calculations Eq. For 
details see text. 
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FIG. 6: Best fit values of Tq for different nuclei (symbols). The dashed line corresponds to Eq. 1361 
used for Z > 15. 
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FIG. 7: Caloric curves are shown for calculations of the system of A=168 and Z=75 at fixed breakup 
density and multiplicity-dependent density. The dotted lines are calculated from the SMM85. The 
dashed lines result when empirical binding energies are taken into account. The solid lines are 
obtained from the improved model, ISMM, with empirical modifications of both binding energies 
and free energies. 
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FIG. 8: Pressure curves due to kinetic motion and Coulomb interaction (see Eq. 140(1 are plotted 
for the system of A=168 and Z=75 at fixed breakup density and multiplicity-dependent density. 
The dotted lines are calculated from the SMM85 while the ISMM presents the solid lines. 
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FIG. 9: Average breakup multiplicities are shown for the system of A=168 and Z=75 at fixed 
breakup density and multiplicity-dependent density. The dotted lines are calculated from the 
SMM85 while the ISMM presents the solid lines. 
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FIG. 10: Dependences of temperature on excitation energy and breakup density are shown for the 
system of A=168 and Z=75. Calculations as function of excitation energy at fixed density of 1/6 
normal density are shown as solid circles in the left panel. Calculations as function of density at 
fixed excitation energy are shown as open squares in the right panel. 
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FIG. 12: Final mass and charge distributions after applying MSU-DECAY, the empirical secondary 
decay procedm'e discussed in Sect. 11111 The dashed lines are calculated from the primary results 
of the SMM85 while the sohd hues are from ISMM. 
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FIG. 14: Primary isotopic distributions for Be, C, O and Ne nuclei. The dashed lines correspond 
to the calculations of the SMM85 while the solid lines represent the results of ISMM. 
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FIG. 15: Isotopic distributions are shown for isotopes from Li to O. Experimental data are shown 
as the solid circles. The dashed lines denote calculations from the SMM code used in ref-fj] and 
the solid lines are the final distributions obtained using the present ISMM model, which contains 
an empirical secondary decay procedure. 



46 



10 
7 

5 



> 2 

llO 
7 
5 

3 
2 



" " '^t^'^fi^) ISMM 
P T('^'^^O) 


T(^-^He) ' 7 : 
■ / \ ' ^ 


' ■ ' ' 1 ■ ' ' ' 1 ' ' 


io_viiVi — prs ec|Uii 
^ ^ 


• 




■ .... 1 .... 1 . . 




. . . . 1 . . . . 1 . . ■ 



10 



15 



5 
Ai 



10 



15 



FIG. 16: Isotopic temperatures extracted from three types of thermometers (see tabled). Experi- 
mental data are shown as the symbols. The lines are the calculations by the ISMM. For reference, 
the primary temperature of 5 MeV calculated from the ISMM is shown as the horizontal dotted 
lines. For details see text. 
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